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ABSTRACT 

It is well-known that the equations for a simple fluid can be cast into what 
is called their Lagrange formulation. We introduce a notion of a generalized 
Lagrange formulation, which is applicable to a wide variety of systems of 
partial differential equations. These include numerous systems of physical 
interest, in particular, those for various material media in general relativity. 
There is proved a key theorem, to the effect that, if the original (Euler) system 
admits an initial-value formulation, then so does its generalized Lagrange 
formulation. 
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1 Introduction 



Consider a simple perfect fluid in general relativity. That is, fix a space- 
time — a 4-dimensional manifold M with metric gab of Lorentz signature 
(—,+,+,+). The fluid is described thereon by two fields, a unit timelike 
vector field (which is interpreted as the velocity field of the fluid), and a 
scalar field p (which is interpreted as its mass density). These fields must 
satisfy the fluid equations, 

(p + p)n™V^M'^ = -{g''^ + u'^u^V^p, (1) 

V^(pm"^) = -pV^M™. (2) 

Here p is specified as some fixed function of p, the function of state. 

This treatment is usually called the Euler formulation of a fluid. Its 
characteristic feature is that the fluid is described by means of fields on 
space-time. That is, the "independent variable" in this formulation — the 
thing the fields are functions of — is the event of space-time. There is an 
alternative treatment of a fluid, called the Lagrange formulation, in which 
we "move with the fluid, rather than remain fixed in space-time". In other 
words, the independent variable in this formulation is the fluid-element, and 
so the fluid is described by fields that are functions on the manifold of fluid- 
element sQ 

Each of these two formulations has its advantages. The Euler formulation 
is less tightly tied down to the fluid itself, and so is usually more convenient 
when other systems — which would naturally be described with reference 
to space-time — are involved. In particular, the Euler formulation is nor- 
mally used for a fluid in interaction with other fields, as, for example, in 
the Einstein-fluid system. The Lagrange formulation, by contrast, tends to 
be more convenient when one wishes to identify and follow individual fluid 
elements. For example, the Lagrange formulation might be used to describe 
a fluid object with a boundary. The boundary, in this formulation, would be 
fixed once and for all at the beginning (by designating those fluid-elements 
that constitute the boundary) as part of the kinematical structure. In the 
Euler formulation of such an object, by contrast, the boundary would be 
"dynamical" . 

^See, for example, for the Euler and Lagrange formulations of non-relativistic perfect 
fluids, and Appendix A of |2] for the Euler formulation of a relativistic perfect fluid. 
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How are the Euler and Lagrange formulations related to each other? Cer- 
tainly, the two are physically equivalent, i.e., they represent mere mathemat- 
ical reformulations of the same physics. That is, all physical predictions will 
be the same, no matter which formulation is used; and, at least in principle, 
either formulation could be used to solve any given problem. Indeed, one 
might be tempted to go further than this, to view them as related by a mere 
coordinate transformation on the manifold of independent variables. But 
such a viewpoint would be misleading, for the "coordinate transformation" 
between the two sets of variables involves the dynamics of the system. Thus, 
for example, from the standpoint of the Euler formulation the Lagrange for- 
mulation represents a curious mixing of kinematics with dynamics. 

These mathematical differences in fact go even deeper. It is well-known 
that the equations for a perfect fluid in the Euler formulation, Eqns. (|l])-(@), 
have a well-posed initial-value formulation^]. But the corresponding equa- 
tions in the Lagrange formulation — at least, those obtained directly, by 
simply "transforming" the Euler equations — do not^. However, it has been 
shown by Friedrich, in |Q, that, at least for a certain fluid system in gen- 
eral relativity, there can be introduced a Lagrange formulation having also 
an initial- value formulation. It is necessary, in Friedrich's treatment, to in- 
troduce a substantial number of additional fields (including a frame-field) 
together with additional equations on those fields. What is not so trans- 
parent, however, is the mechanism behind this treatment. Precisely what 
features of these fluid systems are needed for such a deterministic Lagrange 
formulation? 

The purpose of this paper is to introduce and explore a certain, broad, 
geometrical setting for the Lagrange formulation of systems of partial differ- 
ential equations. 

In Sect 2, we introduce that setting. Our framework is systems of partial 
differential equations that are first-order and quasi-linear (i.e., involving only 
first derivatives of the fields, and those only linearly) — a framework that 
includes virtually every partial differential equation in physics. Given any 
such system — provided only that it has among its fields a distinguished vec- 
tor field — we write out a new system, its "Lagrange formulation" . The key 

^For the case of the Einstein-Euler system, for example, see Sect. 4.2 of [|, and 
references therein. 

■^In fact, some care must be taken, in the Lagrange formulation, even to say what 
"initial-value formulation" means, in light of the fact that the independent variables are 
not the usual space-time events, through which evolution normally proceeds. 
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idea of this scheme is what one might expect: Include, among the dynami- 
cal variables of the new system, what were the independent variables of the 
original system. It turns out that, in order to execute this scheme, it is nor- 
mally necessary to introduce additional dynamical variables and equations. 
We give a general scheme for choosing these variables. The key result of this 
section is the following: Given any system of partial differential equations 
having a distinguished vector field as above, and also having an initial- value 
formulation, then a certain version of its Lagrange formulation also has an 
initial-value formulation. 

In Sect. 3, we give some examples of this scheme. We apply the present 
scheme not only to ordinary fluids, but also to various other types of material 
systems, including dissipative fluids and elastic solids. This scheme is also 
applicable when such material systems are undergoing interaction, e.g., when 
they are coupled to an electromagnetic or gravitational field. Finally, we show 
in Sect 3 how Friedrich's original system fits within the present framework. 

A number of related mathematical issues are discussed in the appendices. 
In Appendix A, we describe a general procedure for modifying any system 
of partial differential equations by "taking derivatives" of the fields of that 
system. This procedure, it turns out, is crucial for casting systems into a 
form in which our Lagrange formulation can be applied. In Appendix B, we 
review a few facts about the initial-value formulation of systems of partial 
differential equations. (For a more detailed treatment, see, for example, 0.) 

2 Lagrange Formulation 

Fix a first-order, quasilinear system of partial differential equations. That 
is, let there be given a fibre bundle, consisting of some base manifold M, 
some bundle manifold B, and some smooth projection mapping B A M. 
Typically, M will be the 4-dimensional manifold of space-time events (but it 
could be any smooth manifold). By the fibre over a point x of M, we mean 
the set of all points y of B such that 7r{y) = x. Think of the fibre over x E M 
as "the set of possible field- values at x" . Then B is interpreted as the set of 
"all possible choices of field- values at all points of M" , and vr as the mapping 
that assigns, to each such choice, the underlying point of M. Thus, point y 
of B could be written as y = {x, 0), with x e M and in the fibre over x. 
The action of the projection mapping would then be given by 7r(a;, 0) = x. 
Typically, the fibre over a point x E M will be some collection of tensors, 



4 



with given index structure (possibly subject to various algebraic conditions), 
at X, whence B will be a manifold of all such tensor-collections at all points 
of M. In this case, B is called a tensor bundle. However, B could in general 
be any smooth manifold, subject only to the local-product condition in the 
definition of a fibre bundle 0. 

By a cross-section of such a bundle we mean a smooth mapping M^B 
such that 77 o is the identity map on M. In other words, a cross-section 
assigns, to each point x of M, a point of the fibre over x; i.e., it assigns a 
"field-value" at each point of M. In the case of a tensor bundle, a cross- 
section is simply a certain collection of smooth tensor fields on M. Our 
partial differential equation will be an equation on this map, linear in its 
first derivative. In order to write out this equation, we introduce two smooth 
fields, k'^°'a and j"^, on B. Since these are fields on B, they depend on point 
y = (x, 0) of B, i.e., they depend on a choice of "point x of the base manifold, 
as well as field- value at that point". The index "a" on k^^'a is a tensor 
index in B at the point, y E B, at which this field is evaluated; the index "a" 
is a tensor index in M at the corresponding point, 7r(?/), of the base manifold. 
The index "A", on both fc^^^Q, and j"^, lies in some new vector space (which 
will turn out, shortly, to be the vector space of equations). Finally, our 
partial differential equation, on a cross-section 0, is 

fc^^(V0)„'^=/. (3) 

This equation is to be imposed at each point x G M, with the fields k and 
j evaluated at 0(a;) G B, i.e., on the cross-section. Here, (V0)a" denotes 
the derivative of the map (i.e., a map from tangent vectors in M at a; to 
tangent vectors in B at 0(a;)). The index "A" in Eqn. (|^) is free, i.e., Eqn. 
(0) represents a number of scalar equations equal to the dimension of the 
vector space in which "A" lies. 

Here is an example. Fix a 4-dimensional manifold M, together with a 
Lorentz-signature metric gab on this M. Let B be the 8-manifold consisting 
of triples, {x,u°',p), where x is a point of M, m° is a unit timelike vector at 
X, and p is a number. Let 'n'{x,u"',p) = x. This is a fibre bundle; in fact, a 
tensor bundle. The fibre over a point x G M consists of (u",p), a vector at 

^Recall that this condition requires, essentially, that, locally in M, B can be written 
as a product, M x F, of Af with some other fixed manifold F, in such a way that the 
projection mapping tt becomes the projection to the 7\f-factor in this product. This 
condition guarantees, e.g., that, locally, all the fibres of the bundle are diffeomorphic with 
this fixed manifold F, and so with each other. 
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X together with a number. A cross-section of this bundle is represented by 
smooth fields, and p, on M. Let the equations, on such a cross-section, 
be where p{p) is some given, fixed function of one variable, and Va is 

the derivative operator defined by the space-time metric Qab- This is a first- 
order, quasihnear system of partial differential equations, i.e., the equations 
are linear in the first derivatives of the fields. The vector space of equations, 
in this example, has dimension four. This system, of course, describes a 
simple perfect fiuid in general relativity. 

We shall now introduce a technique that transforms a given first-order, 
quasilinear system of partial differential equations — provided that system 
lies within a certain class — into a new first-order, quasilinear system of 
partial differential equations. This new system will be called the Lagrange 
formulation of the original. While the new system will differ in many respects 
from the original one — e.g., it will have a different base manifold, a different 
bundle manifold, and a different number of equations — the two will be 
intimately related to each other. In particular, it will turn out that there is 
a natural correspondence between the solutions of the original system and 
those of its Lagrange formulation. 

In order to apply this technique to a given system of equations, it is 
necessary that that system satisfy the following condition: Among the various 
fields of the system there must be distinguished one consisting of a nowhere- 
vanishing vector field on the base manifold M. This condition means, then, 
that the fields of our system take the form {u°',ip), where represents the 
nowhere- vanishing vector field on M, and if represents "the rest of the fields" . 
Thus, given a system that has, among its various fields, no vector field at 
all, then we shall be unable to write out any Lagrange formulation for it; 
and if it has several vector fields, then we must, at this stage, distinguish a 
particular one. We shall denote hj B ^ M the bundle in which the rest of 
the fields, the ip, lie, and use Greek indices for tensors in the manifold B. 
Note that these are different from the Greek indices, e.g., in Eqn. (^, for 
tensors in the manifold B. The equation for our system may now be written 
as 

k'^WaU' + k"^\{V^)a'' = J^, (4) 

where k'^°'b, k"^°'a, and are all functions of u"",^, and point of M. In 
Eqn. (^), the Va in the first term can be any derivative operator on M; 
and the form of j"^ depends, of course, on what operator has been chosen. 
We could, for example, simply fix, once and for all, some derivative operator 
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Va, and use it to write Eqn. (^. Should it happen that the manifold M 
comes equipped with a kinematical metric (i.e., one not included among the 
physical fields then it is often convenient to use its derivative operator 
in Eqn. (^. This possibility is available, e.g., for systems representing fiuids 
in special relativity, or in general relativity with fixed background metric. In 
fact, we could even choose the derivative operator Va in Eqn. (|) to depend 
on the fields (u", if) themselves, provided only that its dependence on these 
fields is algebraic, rather than through their derivatives. We now obtain the 
Lagrange formulation of this system. 

For the base manifold of the Lagrange formulation, we choose any mani- 
fold M having the same dimension as M. Tensors over this M will be denoted 
by lower-case Latin indices with hats. We also fix, once and for all on this 
manifold M, a nowhere- vanishing vector field, u"'. This is a purely kine- 
matical object, i.e., it is fixed right at the beginning, and will not be subject 
to any dynamical equations. 

We next specify the bundle manifold, B, of the Lagrange formulation. 
Fix a point, x, of the base manifold M. Let the fibre over this point consist 
of a triple, {x, ip, Kq,^), where i) x is a point of M, the base manifold of the 
original system, ii) ip is a point of the fibre over x in B, the bundle manifold 
for the original system, and iii) kq,^ is an invertible tensor, where the index 
"a" refers to tensors in M at the point x G M and the index "6" refers to 
tensors in M at the point x G M. A more detailed discussion of these three 
objects follows. 

i) The points (x) of the base manifold M of the original system become, 
in its Lagrange formulation, field-values. In the case of a simple perfect fiuid, 
for example, each point of the original base manifold M represents an event 
of space-time; while each point of the new base manifold M represents "a 
particular fiuid-element at a particular moment of its life". Thus, in the 
Lagrange formulation of such a fiuid, x will be a field over x, a field that 
specifies "which event in space-time that particular fiuid-element occupies at 
that particular moment" . 

ii) The field-values, the (f, of the original system become field- values also 
in its Lagrange formulation. But there is one important change: What were 
fields over M in the original system become, in its Lagrange formulation, 
fields over M. Thus, were the fields collected in ip all tensor fields on M, 
then the corresponding fields in the Lagrange formulation would depend on 
point X of M, but would continue to be tensors in the tangent space at the 
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point X of M0. In the case of a simple perfect fluid, this step amounts, 
physically, to "attaching the density p to the fluid element, rather than to 
the point of space-time." 

iii) There is introduced a new object, Ka^, an invertible two-point tensor, 
with one index at a; G M, the other at x G M. Nothing analogous was 
present in the original system. Denote the inverse of /ta^ by k;,'^, so we have 
i^a^b" = ^a" and Kft^Ka'^ = h^- The role of this tensor Ka'' is, as we shall see, 
to preserve the first-order character of the final system of equations. Note 
that the dynamical field in the original system has disappeared entirely: 
There is no analog of it as a dynamical field in the Lagrange formulation. 

Next note that the pair (x, where x is a point of M and <y9 is a point 
of the fibre in B over x, is precisely the same thing as a point of the bundle 
manifold B. Call that point (for later convenience) so we have = 
(x, v?) G B. Then we may recover the point x of the original base-manifold 
from the point & B using the projection vr: We have x = ^{(p). Thus, our 
construction of the bundle manifold B for the Lagrange formulation could 
have been stated as follows: The fibre over point x G M consists of a pair, 
[0, Ka'^), where is a point of the manifold B, and Ka'' is an invertible tensor 
with one index at x G M, the other at vr(<^) G M. 

We have now completed the specification of the fibre bundle in which the 
Lagrange formulation of our system will be written. The base manifold, M, is 
some new manifold, of the same dimension as M, while the bundle manifold 
B is such that the fibre over x G M consists of a pair, {(f, k^), where & B, 
and Ka^ is a certain 2-point tensor. A cross-section of this bundle, then, is a 
smooth map (a map we also denote by 0) that assigns, to each point x G M, 
a point of B together with a suitable tensor kq,''. On such a cross-section, 
we now impose the following equations: 

(V(7^o<^))a'' = /^a^ (5) 

V[,{Ka]') = ha\ (6) 

k'^^f^J'^iina'u') + k"^\Ka\V0)t^ = j^. (7) 

These are the equations of the Lagrange formulation. In Eqn. (|^), the 
combination vr o is a map from M to M, for goes from M to B, and 

^Note that we can, in this case, convert these to ordinary tensors on the manifold M 
by using and its inverse. This, a mere "coordinate transformation" on the fibres, 
changes nothing, in particular, not the final partial differential equations of the Lagrange 
formulation. 
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vr from B down to M. Eqn. (^) asserts that the derivative of this map 
is precisely the tensor k^''. Thus, this equation provides the geometrical 
meaning of the field kq,^. Note that invertibility of kq.^ in Eqn. (^) implies 
that the map tt o cp from M to M is a local diffeomorphism between these 
two manifolds. It was to achieve this feature that we originally choose M to 
have the same dimension as M. Eqn. is merely the curlj^ of Eqn. (^. 
Any derivative[] may be used on the left in Eqn. (^), but the exact form of 
the function fca^ (of {(p, Ka^)) that appears on the right will depend on which 
derivative was chosen. This situation is analogous to that of Eqn. (^. Eqn. 
(|^) is the translation of the equation of the original system, (||), to our new 
system. Here, everywhere in the fields fca^, k'^""},, k"^"'a, and there is to be 
substituted the combination "Ka^'w^" for "m*"'; and "<^" for 'V"- In Eqn. (0), 
this "replacement" takes place even inside the derivative. Note that the field 
of the original system has now disappeared entirely, having been replaced 
by the image of the kinematical field under the mapping ir o (p. 

Thus, beginning with any first-order, quasilinear system of partial differ- 
ential equations of the form (^, we obtain a new system of equations, its 
Lagrange formulation, of the form ©-(l^). The Lagrange formulation has a 
completely new base space, but fields and equations that echo those of the 
original system. 

We now claim: Every solution of the Lagrange formulation gives rise, at 
least locally, to a solution of the original system. Indeed, let {(p, k^^) be fields 
satisfying (^)-(0). Then, as we have seen, it o ip is a. local diffeomorphism 
between M and M. We now introduce the following two fields on M: u'' = 
{'V{7!'o0))iJ'u"', and ip = 0o (iio0)~^. That is, we let and ip be the images 
of u'' and 0, respectively, under the diffeomorphism tt o 0. Then these fields, 
(m^, ip), on M satisfy the system (^), as is immediate from Eqns. (|^),(0). We 
next claim that the converse also holds: Every solution of the original system 
gives rise, at least locally, to a solution of its Lagrange formulation. Indeed, 
let {u'',ip) be fields satisfying (^. Choose any manifold M with the same 
dimension as that of M, and any nowhere- vanishing vector field thereon. 

^For convenience, we shall always include within our system all first-order equations 
on the fields of the system, even those that arise from differentiating other equations of 
the system. 

''These derivatives may be characterized in the following manner. Consider the bundle 
with base space M and fibre over x € M consisting of a pair, {x, Kq.^), where x € M and 
Ka is a tensor with indices at x and x. Then a choice of connection in this bundle gives 
rise to an operator Va for use in the left side of Eqn. (||). 
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Now let 0he a diffeomorphism between M and the cross-section, (^[M], such 
that (vr o 0) sends u to u; and then define by Eqn. (||). Then these fields 
{ip,Ka^) on M will satisfy Eqns. (the first two by construction, the 

last by Eqn. (D). 

Thus, the original system and its Lagrange formulation are identical as 
to solutions. But the two systems are quite different as to form. Their base 
manifolds, M and M, although of the same dimension, differ in their geom- 
etry. The manifold M must be endowed with a fixed, kinematical "velocity 
field" , m", while M has no such kinematical field. On the other hand, various 
kinematical fields that might have been specified over M (such as a Lorentz 
metric) yield no analogous kinematical field^ on M. Furthermore, the fields 
of the two systems differ in several respects. Beginning with the fields of the 
original system, we must delete the dynamical field m", while adding "point 
of M" as well as the invertible tensor Ka^, to obtain the fields of the Lagrange 
formulation. Finally, the equations for the two systems differ in that, for the 
Lagrange formulation, there must be introduced one new equation, (^, on 
the derivative of the "point of M", as well as is the curl, (|^), of this new 
equation. 

What we have described above is precisely what is usually done in writing 
down the Lagrange formulation for a fluid. For example, consider again 
the simple perfect fluid, with flelds (m'^,p) on M and equations (|I])-(^. Its 
Lagrange formulation consists of flelds^ {x,K,f^"',p) on M, with equations 
consisting of (^, (^, and 

(p + p{p))u'V,iK^,^un + (^7"'" + M Vm"/€^"^)/€„^V^p(p) = 0, (8) 

w'VgP + (p + p(p))/€a'Vs(/€™'^M") = 0. (9) 

We now return to the general case. It turns out that the procedure given 
above — starting with a system and ending with its Lagrange formulation 

®A Lorentz metric on M, for example, becomes, on AI, an algebraic fimction of the 
fields (namely, just of x) of the Lagrange formulation. 

^In the notation of (|5|)-(^, we have (p = {x, p), and ip — {x, p). 

^"There is an unfortunate complication here, involving the normalization condition, 
vf^u^Qab = — 1, on u". It is awkward simply to carry this condition through the La- 
grange formulation. But there are several other ways — none very elegant — to deal with 
it. Perhaps the simplest is to rewrite the fluid equations from the outset (by inserting, 
strategically, factors of {u°'u^gab)) in such a way that, while retaining their initial-value 
formulation, they no longer require this normalization condition. Then take the Lagrange 
formulation of these new equations. 
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— suffers from a serious difficulty. In general, the equations of the Lagrange 
formulation, will fail to have an initial- value formulation, even if the 

original system, (^, did have such a formulation. For example, the system 
(H)"®) di)"® has no initial- value formulation, although the system (|l|)-(0) 
of course does. But it turns out that this difficulty does not arise — i.e., 
the Lagrange formulation does inherit an initial-value formulation from the 
original system — provided the original system satisfies the following condi- 
tion: There can be derived from Eqn. (^) an expression for the derivative of 
the vector field m", without contractions, back in terms of the various fields 
of the system. In other words, it must be possible to cast Eqn. (|[) into the 
form 

V.m" = Wa\ (10) 
fc"^\(V^)a"=/^, (11) 

where wj', k"^"'a and j'^ are functions of (x, u"", ip), i.e., are functions of the 
point of B and the vector u"'. In Eqn. (p^, Va can, again, be any derivative 
operator on the manifold M; and the form of wj' depends, of course, on what 
operator has been chosen. Note that, once we have derived from Eqn. (^) an 
equation of the form (0), then it is easy to cast the equations that remain 
into the form (|11|): Simply use Eqn. (p!0| ) to remove all M-derivatives from 
Eqn. (H). Indeed, we have j'^ = — k'^"-},Wa. 

The equations for systems of physical interest typically do not take the 
form of Eqns. ([TO|)-([lID, i.e., they do not express the derivative of u"" in terms 
of the other fields. For example, Eqns. (|l])-(@) do not have this form. But it 
turns out that there is a simple, general procedure by which any first-order, 
quasilinear system of partial differential equations having a preferred vec- 
tor field u"" can be recast so as to take the form (|T0|)-(|11]). This procedure, 
called taking the derivative system, is spelled out in Appendix A. It consists 
of modifying the original system by introducing additional fields, which rep- 
resent the derivatives of the original fields, as well as additional equations 
on those fields. The result of taking the derivative system is to produce a 
new system of partial differential equations, having, in an appropriate sense, 
identical solutions to the original. Applied to a system in which a preferred 
vector field has been distinguished, it produces a system in which Vqm'' 
is expressed back in terms of the fields of the system. Furthermore, applied 
to any system having an initial-value formulation, the derivative system also 
has an initial-value formulation. 
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As an example of this procedure, we return to the system, for 
a simple perfect fluid in general relativity. For the distinguished nowhere- 
vanishing vector field in this case, we choose, of course, the velocity field 
u"' of the fluid. The result of taking the derivative system of this system is 
the following. The fields consist of (u", p, wj', Va), where is a unit timelike 
vector field, p a positive scalar field, Wa'^ a tensor field satisfying QabU^'wJ' = 0, 
and Va a vector field, all subject to the algebraic conditions 

(p + p)u"'w^^ + {g'^"' + {dpldp)v^ = 0, (12) 

u'^Vm + {p + p)w^'^ = Q. (13) 

On these fields is imposed the following system of first-order, quasilinear 
partial differential equations 

Van" = Wa\ (14) 

V^aW^'' = Rabm'u"', (15) 
VaP = Va, (16) 
"^[aVb] = 0. (17) 

Note what has happened here. We have introduced two new fields, Wa'' and 
Va- The "interpretation" of (via (0)) is as the derivative of and of Va 
(via (|1^)) as the derivative of p. The original fluid equations, (|l|)-(|^), have 
been converted into algebraic conditions, ([l^)- ([T3|) , on these new fields. That 
is, the original fluid equations serve merely to define the bundle of fields for 
this new system. Finally, the new system contains two other equations, Eqns. 
( [T5| ) and ([T7|), that are merely the curls of Eqns. (|T3|) and (|16|), respectively. 

In short, our "procedure" has done nothing of substance. But note that, 
starting with a system, (|lj)-(^, which fails to express Vqm" in terms of the 
fields of the system, our procedure produces a new system satisfying, via 
([T^), this condition. Furthermore — and this is perhaps the striking feature 
— the system (p!^)-([T7|) inherits from the original fluid system, (|l|)-(||), its 
initial-value formulation. 

The key result of this section is the following: Consider any system, 
of partial differential equations in which there has been selected a preferred 
vector field u"" . Let i) that system have an initial-value formulation, and 
a) the equations of that system express the derivative of u'^ in terms of the 
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fields of the system (as in l\n])-(^Tl\)). Then the Lagrange formulation of that 
system also admits an initial-value formulation. 

First note that the Lagrange formulation of the system (|l^)-(^) consists 
of Eqns. d^)-®, together with 

Ra'Vci^Ju'') = Wa\ (18) 

k"^\na\V^)C=f^- (19) 

As discussed in Appendix B, in order that a general first-order, quasilinear 
system of partial differential equations have an initial-value formulation it is 
necessary that it satisfy three conditions: i) the system admits a hyperboliza- 
tion; ii) all the constraints of the system are integrable, and iii) the system 
has the correct number of equations relative to the number of its unknowns. 
What these conditions mean is also explained in Appendix B. We check these 
three conditions in turn. 

Let the original system, Eqns. (|10])-(|11|), admit a hyperbolization. Then 
the construction that, applied to Eqns. ([T0| ) - ([TT1) to obtain a bilinear ex- 
pression in 5ip" yields, when applied to Eqns. (|18|)-(|19|), a corresponding 
bilinear expression in Next, contract Eqn. (|^) with u'^ and use Eqn. 



([TSp to obtain an equation expressing u^Wm i^d algebraically in terms of 
the fields. From this there follows immediately an appropriate bilinear ex- 
pression in 5Ka^. Finally, a bilinear expression in 5x arises from Eqn. (^. 
These three bilinear expressions, taken together, represent a hyperbolization 
for the system (§)-(§), (|TB|)-(|T3). 

Every constraint of the original system, (p!oD -([riD, gives rise to a constraint 
of its Lagrange formulation; and, furthermore, if these constraints of the 
original system are integrable, then so are the corresponding constraints of 
the Lagrange formulation^ This assertion is immediate from the fact that 
Eqns. (|18D and (p!9| ) mimic Eqns. (p!0| ) and (|11|), respectively. But, it turns 
out, there are two additional classes of constraints for the system of the 
Lagrange formulation. The first class arises from taking the curl of each side 
of Eqn. (^). These constraints are certainly integrable, and, indeed, the 
corresponding integrability conditions are precisely Eqn. (^. The second 
class of constraints arises from taking the curl of each side of Eqn. (^. These 
constraints are also integrable, and indeed their integrability conditions are 

^^Note in particular that the original system, (|lo|)- (|ll|), always possesses the constraints 
arising from the curl of Eqn. (|lO|). Thus, if the constraints of this system are to be 
integrable, this curl-equation must have been included in the system (p^-(|lTl). 
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identities, simply from the way Eqn. (j^) was obtained. We conclude, thus, 
that a system of the form (|l^)-(0) having all its constraints integrable leads 
to a Lagrange formulation, d^)-®, (|l8])-(0), also having all its constraints 
integrable. 

Finally, in order to check the third condition, we introduce the following 
integers. Denote by n the dimension of the base space M (the number of 
independent variables of the system), by u the dimension of the fibres in the 
bundle B (the number of unknowns represented by y?), by e the dimension 
of the vector space in which the index "A" of Eqn. ([TI|) lies, and by c 
the dimension of the space of vectors of the form WmC^'A-, as d^A runs over 
constraints for Eqn. (|TT|) . Then, for the original system, we have: the number 
of unknowns is given hy Uq = u + n (the term "n" arising from the field 
the number of equations is given hy Co = + e (these terms arising from 
Eqns. (|TUp and (|TTp, respectively); and the number of effective constraints 
is given by Co = n{n — 1) + c (these terms arising from the constraints of 
Eqns. ([To| ) and (|Tll), respectively). For the Lagrange formulation, on the 
other hand, we have: the number of unknowns is given hy = u + n + 
(the term "n" arising from the field "point of M", the term "n^" from the 
field Ka'') ; the number of equations is given by = + n^(n — l)/2 + + e 
(these terms arising from Eqns. (||)-(||), (P^ -(P^, respectively); and the 
dimension of the space of effective constraints is given by cl = n{n — 1) + 
n{n — l){n — 2) /2 + n{n — 1) + c (these terms arising from the constraints of 
Eqns. (§P-(§P, (|TBp-(^), respectively). It is easy to check from these formulae 
that Co — Co = Uo implies cl — cl = ul- In other words, if the original system 
has the appropriate number of equations relative to its number of unknowns, 
then so does its Lagrange formulation. 

Thus, we have shown a system of the form (|lOD-(|lI]) having an initial- 
value formulation gives rise to a Lagrange formulation also with an initial- 
value formulation. 



3 Examples 

In this section, we introduce various examples of physical systems, the partial 
differential equations that describe them, and the Lagrange formulations of 
those partial differential equations. 

One such example, the simple perfect fluid, has been discussed already 
in Sect. 2. The fields, on space-time, M,gab, consist of a unit timelike 
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vector field u"" (interpreted as the fluid velocity) and a positive scalar fleld 
p (interpreted as the mass density); and the equations are (|l|)-(||), where 
p(p) is some flxed function (the function of state), which specifles the type of 
fluid under consideration. This is the Euler formulation. In order to achieve 
a Lagrange formulation for this system, the first step is to modify these 
equations so that the derivative of u"", without contractions, is expressed 
in terms of the other fields. This was achieved by taking the derivative 
system: We introduced two new (tensor) fields, Wa and fa, subject to the 
algebraic conditions ([I^)- ([T3|) . We then imposed on the total set of fields, 
(«", p, Wa'', fa), the partial differential equations (|1^-(|1^. This new system 
([T^)- ([T7|) is, by virtue of Eqn. (0), of the required form, and, in addition. 



it inherits from the original system, (|T])-(^, its initial- value formulation. To 
this system, (|l4D-(p!^), we may therefore apply the methods of Sect. 2 to 
obtain its Lagrange formulation. There result fields (x, p, wj'^ Va, i^h') on M, 
subject to the equations (^-(0). This new system, as demonstrated in Sect. 
2, again has an initial- value formulation. 

There is a natural generalization of this simple perfect-fiuid system to 
a much broader class of fiuids. Fix some smooth manifold 5*, the points of 
which will, shortly, be interpreted as representing "local, internal, states of 
the fiuid". Also fix any space-time, {M,gab)- Let the fields, on this space- 
time, consist of a unit, timelike vector field, m° (again interpreted as the 
velocity field of the fiuid), together with a second field, ip, which is valued 
in S (and which is interpreted as giving the local state of the fiuid at each 
point of space-time). Thus, is a mapping, M ^ S. As an example, the 
simple perfect-fiuid system discussed above is the special case in which S is 
a 1-manifold (whose points are labeled by a coordinate p, whence ip reduces 
to the density field p). That is, our simple perfect fiuid is one whose local 
state is completely characterized by the value of the density. 

We next wish to write equations on these fields. To this end, fix two 
tangent vector fields, and T", and one covector field. Fa, on the manifold 
S, where we have introduced Greek indicesQ to represent tensors in S. The 
physical interpretations of these fields will be given shortly. Let the equations 
for this system be 



u^Vau' + ig"' + nV)(V^)a"F„ = 0, (20) 



^^These are not to be confused with the indices for tensors on the bundle space, used 
extensively in Sect 2. 
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M"(V(/?)a" + V'VaU'' + T° = 0. 



(21) 



The first equation gives the fluid acceleration in terms of the derivative of the 
fluid state. We may interpret the fleld F^, which acts by driving the fluid, 
as an "effective force" . The second equation gives the time rate of change of 
fluid state in terms of that state and the divergence of We may interpret 
the flelds and T", respectively, as giving the rate of change of fluid state 
under small volume- changes of a sample of that fluid, and under allowing a 
sample of that fluid to evolve in time. The simple perfect fluid, for example, 
has Fa = {p + aVi = {p + p)d/dp, and T" = (for these choices 

reproduce Eqns. (|l|)-(H)). Another familiar example is the perfect fluid with 
2-dimensional manifold S of internal states, where the additional degree of 
freedom is represented by a conserved particle-number n. In this case, is 
given by the same expression as above, V" by (p + p) d/dp\n + n d/dn\p, and 
again T"' by 0. A more exotic example is that of a fluid consisting of several 
species of particles, between which chemical reactions can take place as the 
fluid evolves. In this case, we would have dim(S') > 2 (the additional degrees 
of freedom describing the chemical composition of the fluid) and nonzero 
(representing the rate and direction of the chemical reactions). 

When does the system above satisfy the three properties, as discussed in 
Appendix B, for having an initial- value formulation? Two of these properties 
are immediate: Clearly, this system has no constraints, and the dimension 
of its space of equations is the same (namely, dim(S') -|- 3) as the dimension 
of its space of flelds. As for the third condition, this system, it turns out, 
admits a hyperbolization if and only ifQ V^Fa > everywhere on S. Note 
that, in the explicit examples given above, the combination V^Fa is precisely 
the square of the sound speed. 

We now have a system of equations, (p0|)-(|21|), having a preferred vector 

^■^Note that the last two terms on the left in Eqn. ( pl| ) constitute the most general 
expression (involving and ip) quasilincar in the derivative of m". 

^*For "if , suppose that V"Fa > everywhere on S. It follows that there exists a 
positive-definite metric field, gap, on the manifold S such that V^gap — Fp everywhere. 
Choose one, (e.g., the sum of FaFp / {F^V^) and a suitable positive semi-definite tensor 
haf3 that annihilates V^) and consider the bilinear expression 

- (w^w") [gabSu'^S'v!' + gc^pSf'^S'ipf'] - [5u™5V" + S'u"'6(p°'] . (22) 

This bilinear expression indeed arises, as described in Appendix B, from Eqns. (pO[)-(pT|), 
and is indeed positive-definite (for Wm sufficiently close to Um)- So, this bilinear expression 
gives rise to a hyperbolization. The converse is easy. 
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field, m", and, subject only to the inequality V^'Fa > 0, having an initial- 
value formulation. So, we may apply to this system the results of Appendix 
A and Sect. 2. The first step is to take the derivative system (Appendix A). 
The result of this step is to include, in addition to the fields m", (p above, two 
new fields, Wa'' (with UbWa'' = 0) and Ca", subject to the algebraic conditions 
m'^u7„*+ (c/'*^ + uV)Ca"i^a = and uXa" + V"Wa'' + T'' = 0. (These algebraic 
conditions reflect Eqns. (pOD-(^).) The equations on these fields for the 
derivative system are given by 

VaU' = Wa\ (23) 
V[„Wfef = Rabdu\ (24) 

(V^)," = Ca", (25) 
V[aC6]" = 0. (26) 

This system indeed has a preferred vector field, u"; has among its equations 
one (Eqn. (^)) that expresses the derivative of this algebraically in terms 
of the fields; and has an initial-value formulation (by virtue of that for Eqns. 
([20|)-(|21|)). So, we may, as described in Sect 2, take the Lagrange formulation 
of this system. There results a new system of partial differential equations, 
(|^)-(|^), again having an initial- value formulation. 

Even the broad class of generalized fluids above does not include all pos- 
sible types. For example, there exist fluids manifesting dissipative effects, 
such as heat-flow and viscosity. One description of such a fluid in relativity 
([H) [H) 0! [H) proceeds as follows. The fields consist of a unit timelike 
vector field u"" (interpreted as the fluid 4- velocity), two scalar fields, p and 
n (interpreted, respectively, as the fluid mass density and particle-number 
density), a vector field satisfying u^-qa = (interpreted as the heat-flow 
vector), and a symmetric tensor field Tab satisfying u'^Tab = (interpreted 
as the stress tensor). Thus, the space of field- values at each point of M is 
14-dimensional. The equations on these fields consist of i) vanishing of the 
divergence of nu"" (conservation of particle number), ii) vanishing of the di- 
vergence of {p + p)u°'u'' ^pg"'^ + 2u^"'q''^ + r"^* (conservation of stress-energy), 
and iii) a certain system of nine additional equations that, effectively, gov- 
erns the dynamical evolution of and r"^. It turns out that the resulting 
system, consisting of i)-iii), has an initial- value formulation: Specifically, it 
has a hyperbolization and no constraints. Furthermore — and this is perhaps 
surprising — this system of equations can be so chosen that it reduces, in an 
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appropriate limit, to the familiar Navier-Stokes system for a dissipative fluid. 
(The Navier-Stokes dissipation coefficients — the thermal conductivity and 
viscosity — arise from within the nine equations iii).) Here, in any case, is a 
system of equations with a preferred vector field u"" — a system, therefore, 
to which the present methods can be applied. Thus, we take the derivative 
system, as described in Appendix A, and then the Lagrange formulation, as 
described in Sect. 2. There results a Lagrange formulation for a dissipative, 
relativistic fluid. 

There exist still other types of material systems, e.g., some that are not 
fluids at all. Consider, for example, the elastic solid. In one treatment^ 
of such a system in relativity, the fields consist of a unit timelike vector 
field (the material 4- velocity), a positive function p (the mass density of 
the material), and a symmetric tensor field hah satisfying habv!^ = 0. This 
hab represents the geometry of the material as it was "frozen in" at the time 
the material originally solidified: It describes the shape to which the material 
would "like to return" . Thus, the combination hab — {gab+UaUb), the difference 
between this natural geometry and the actual spatial geometry in which 
the material currently finds itself, is interpreted as the strain of the solid 
material. The equations on these fields are Cuhab = (the vanishing of the 
Lie derivative of hab, interpreted as asserting that the material remembers, 
over time, its frozen-in geometry), and Vb{pu"'u'' + t"-^) = 0, (interpreted 
as conservation of stress-energy, whence r**^ is interpreted as the stress of 
the material). Here, 

^ab 

to be given as some fixed function of hab, Qab, 
and m". This is the stress-strain relation. Provided this stress-strain relation 
is chosen appropriately, the final system, it turns out, has an initial-value 
formulation: Specifically, it has a hyperbolization and no constraints. Again, 
we have a system to which the present methods can be applied. There results 
a Lagrange formulation for an elastic solid. 

There are, presumably, a variety of other systems of equations, repre- 
senting "materials" of various sorts, having, among their fields, a preferred 
4- velocity. Examples might include the systems for a plasma, for a supercon- 
ductor, or for a solid (such as ice) that is able to flow. These systems, too, 
will have Lagrange formulations. 

These various material systems may, of course, interact with their en- 
vironment in a variety of ways, e.g., electromagnetically, gravitationally, or 
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For other treatments, as well as the local existence theory for solutions, 



see [p| and ||10|. For a brief summary of this subject, see 
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through contact forces. What impact do such interactions have on their 
Lagrange formulations? 

Consider, as an example, the fluid of Eqns. (pOD-(pTD interacting electro- 
magnetically. This charged-fluid system is described by fields consisting of 
the original fluid variables, and together with an antisymmetric (elec- 
tromagnetic) tensor field Fab- The equations on these fields consist of Eqn. 
([20|), modified by the inclusion of a term on the right of the form fiF^aU"", 
Eqn. ([21|)p|, and Maxwell's equations, 

V'Fab = aua, (27) 

V[,Ffe,] = 0. (28) 

Here, the /i in the first equation and the a in Eqn. (^) must be given as fixed 
fields on the manifold S of fluid states. The field a describes how the fluid 
drives the electromagnetic field, and so is interpreted as the charge density. 
We require that it satisfy charge conservation: ^"VaO" = cr, T^VqCt = 0. 
The field /i, which describes how the electromagnetic field drives the fluid, 
might be called the specific charge density. (For a normal fluid, a and fi are 
in ratio (p + p)-) Here, in any case, is a list of fields, together with a system 
of equations on those fields. This system has an initial-value formulation, 
which it inherits from the separate initial-value formulations for the original 
fluid system ((^D|)-(^T])) and for Maxwell's equations. We wish to take the 
Lagrange formulation for this system. Since the system does not express 
the derivative of m° in terms of the other fields, the first step is to take 
the derivative system. But note that, in taking the derivative system, it is 
necessary to introduce, not only the new fields wj' and Ca" that represent 
(via Eqns. (|23| ) and (pSj), respectively) the derivatives of and (p, but also 
the field (ate that represents (via Eqn. (^If)) the derivative of Fab- One might 
have hoped that it would be possible, exploiting somehow the fact that our 
system of equations splits naturally into "fluid equations" and "Maxwell-field 
equations" , to avoid introducing the additional field (abc- Unfortunately, this 
seems not to be the case. This issue is discussed briefly in Appendix A. In 
any case, this derivative system has the appropriate form (a preferred vector 
field m", whose derivative is expressed in terms of the fields of the system), 
and an initial-value formulation (which it inherits from that of the original 

^^There could also be included on the right side of this equation terms algebraic in 
the electromagnetic and other fields. Such terms would represent, e.g., an effect of the 
electromagnetic field on the the rates of chemical reactions. 
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coupled system). So, we may apply the methods of Sect. 2. Thus, there is a 
Lagrange formulation for a charged fluid, but it requires the introduction of 
a further field (abc, representing the derivative of the Maxwell field. 

In a similar way, we may write down the Lagrange formulation for a 
charged dissipative fluid, a charged elastic solid, etc. In each of these cases, 
it is necessary to introduce the auxiliary field (abc- 

The situation for gravitational interactions is similar. Consider, again, the 
fluid of (pOD-(pID, now interacting gravitationally. The interacting system 
is described by fields consisting of the original fluid variables, u"" and ip, 
together with the variables for gravitation: a Lorentz-signature metric gab, 
and a derivative operator, Va- The equations of this system consist of Eqns. 
([2(]|)-(PT|)["1, the equation VaQbc = 0, and Einstein's equation. 

Gab = Tab, (29) 

where Gab is the Einstein tensor. Here, Tab is some fixed symmetric tensor 
function of Qab and the fluid variables (which we interpret as the stress- 
energy tensor of the fluid). It plays a role analogous to that of the functions 
H and a for electromagnetic interactions. We demand of this tensor function 
that, as a consequence of Eqns. (^)-(|2l|), it be conserved[3. This system 
of equations does not have an initial-value formulation, in the sense we are 
using this term. But this is merely a consequence of the fact that our sense of 
this term is overly restrictive, in that it does not tolerate the diffeomorphism 
freedom characteristic of all systems in general relativity. In a physical sense, 
i.e., once the diffeomorphism freedom has been treated properly, the fluid- 
Einstein system does, of course, have an initial-value formulation. Now take 
the derivative system of this system. Note that in doing so we must, as 
in the electromagnetic case, include also fields to represent the derivatives 

^^Note that there are no expressions, algebraic in the gravitational fields, that could 
be introduced on the right in these equations. This is a reflection of "the equivalence 
principle" . 

^^The most general candidate for such a stress-energy (i.e., the most general algebraic 
function of our fields, having the correct index-structure) is given by Tat = {p + p)uaUb + 
Pdab, where p,p are some functions on the manifold S of fluid states. When does there 
exist such a T°-^ that, in addition, is conserved, VbT"^^ = 0, by virtue of the field equations 
(pO|)-(^)? It is not difficult to check that (assuming V^F^ > 0; and demanding p + p > 0) 
a necessary and sufficient condition is that the fields Fa, V", and T" on S satisfy the 
following three equations: i^[aV^F^] = 0, T^Ka = 0, and \7[a{Ki3j + Fp]) = 0, where we 
have set = (2y'3V[^F„] + Fa)/{V^F^). 
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of the gravitational fieldsQ. Take the Lagrange formulation of the result. 
The resulting system, again, will not have an initial-value formulation in our 
restrictive sense, but it will have such a formulation when the diffeomorphism- 
freedom is properly taken into account. We conclude, then, that there does 
exist a Lagrange formulation for a gravitating fluid, but that it requires that 
we introduce further fields to represent the derivatives of the gravitational 
fields. 

In a similar way, we may write down the Lagrange formulation for a 
gravitating dissipative fluid, a gravitating elastic solid, etc. In each case, it is 
necessary to introduce fields representing the derivatives of the gravitational 
fields; and in each case the Lagrange formulation retains the initial-value 
formulation of the original system. 

A similar treatment is available for systems consisting of two or more 
different materials in interaction. In these cases, there will be two or more 4- 
velocity fields present, and we shall have to select one to be that with respect 
to which the Lagrange formulation is taken. 

The treatment of systems in which several interactions are turned on 
simultaneously, e.g., the charged gravitating fluid, is similar. 

Finally, we briefly characterize, within the present framework, Friedrich's 
Q original example of a relativistic Lagrange formulation. Begin with the 
system for a gravitating fluid, as described above, for the case in which 
the fluid has a 2-dimensional manifold S of local states, i.e., that in which 
T" = and V^" = {p + p)d/dp\n + nd / dn\p. For this system, first take 
the derivative system, and then the Lagrange formulation. The result of 
this process — after three, essentially cosmetic, further modifications — is 
precisely Friedrich's original example. The three further modifications are 
the following. 

1. Introduce, already in the original Einstein- fluid system, before taking 
the derivative system, a 3-dimensional space of additional variables, consist- 
ing of three unit vector fields, x'*, y", and z"", that are required to be or- 
thogonal to each other and to the 4-velocity u"". On these fields, impose the 
equations that they be Fermi-transported by u"". The introduction of these 
fields with these equations does not interfere with the initial-value formula- 

^^In the resulting system, there will initially be two versions of "the derivative of the 
metric gab \ one being the original derivative operator Va, and the other arising (via gab) 
through passage to the derivative system. These two versions are then be set equal to each 
other, via Eqn. A similar phenomenon occurs, e.g., on taking the derivative system 

of the Klein-Gordon system. 
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tion. These fields, which have no direct physical significance, are introduced 
to facilitate the writing of various equations. 

2. After taking the derivative system, but before passing to the Lagrange 
formulation, suppress half of the field Ca", which represents the derivative of 
the fiuid state0. While such suppression of variables will in general destroy 
the initial-value formulation for a system, it turns out that, in this particular 
instance, it does not. Thus, the essential effect of this modification is to 
reduce by four the number of independent variables. 

3. Write the final equations, after passing to the Lagrange formulation, 
not in terms of the specific fields listed above, but rather in terms of others 
that are algebraic functions of these. This choice of variables — choice of 
"coordinates" on the bundle space — is, of course, a matter of convenience. 

4 Conclusion 

We have introduced a scheme that takes a first-order, quasilinear system of 
partial differential equations and produces from it a new first-order, quasi- 
linear system, its "Lagrange formulation". The key requirement, on a given 
system of equations, in order that this scheme be applicable to it is that that 
system have, among its fields, some nowhere-vanishing vector field. Why 
this special role of a vector field? Could, for example, a similar scheme be 
developed based on some other geometrical object (s)? It turns out that there 
are two special features of vector fields that we used in the construction of 
the Lagrange formulation. 

First, nowhere vanishing vector fields on manifolds are locally homoge- 
neous. This means the following. Let there be given any manifold M, any 
nowhere-vanishing vector field u"" thereon, and any point x G M; and, sim- 
ilarly, some other manifold M (of the same dimension), vector field and 
point X G M . Then there always exists a diffeomorphism between neighbor- 
hoods of X and x that sends u"" to u". In other words, nowhere vanishing 
vector fields are "locally all the same": They carry no local structure. We 
used this fact in Sect. 2 in order to replace u"" on M by some kinematical 

^°This "suppression" proceeds, in more detail, as follows. Choose on the 2-nianifold S, 
a function s (which is interpreted in Q as the entropy per particle) satisfying F"VqS — 0. 
Now delete the field Ca" everywhere, by replacing the component Ca^VaS of Ca" by some 
new field fa, and the remaining components of Co" by (V(/3)a"- To the resulting system 
add those further equations that are required for integrability of the constraints. 
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field m'^ on M. 

Second, by virtue of the appearance of the vector field u"" on the left in 
Eqn. ([T8|), the system ffl), (|T8|) for the two-point tensor admits a hyper- 



bolization. We used this fact in Sect. 2 in order to achieve a hyperbolization, 
and consequently an initial- value formulation, for the entire system (H)-®, 

It appears that, given any other geometrical structure manifesting these 
two features, then there could be developed a "Lagrange formulation" based 
on it. It is only necessary to make three key modifications in Sect 2 (all 
involving replacing the vector field by the totality of fields in the new geo- 
metrical structure): i) Replace Eqn. (10) by equations for the derivatives of 
all the fields of the geometrical structure; ii) endow the base manifold M of 
the Lagrange formulation with kinematical fields consisting of all the fields 
of the geometrical structure; and iii) replace Eqn. ( [T^ ) by the corresponding 
equation involving all the fields of the geometrical structure. Unfortunately, 
it is not so easy to find geometrical structures having the two features de- 
scribed above, in part because they are somewhat in opposition to each other: 
The first feature, local homogeneity, prefers fewer fields, relatively devoid of 
structure; while the second feature, hyperbolicity of (^, (p!8|), prefers many 
fields, of rich structure. 

There are a variety of geometrical structures that are locally homoge- 
neous. Examples include: two commuting, pointwise independent vector 
fields; a nowhere- vanishing, curl-free 1-form; a symplectic structure; a flat, 
Lorentz-signature metric. Examples of geometrical structures that yield a 
hyperbolization for (P), ([TsD are somewhat less plentiful. One simple class 
consists of those in which the geometrical structure is comprised of a nowhere- 
vanishing vector field m*^, together with any additional fields of whatever type. 
For structures in this class, a hyperbolization for (^, (^) (suitably gener- 
alized) is guaranteed already by the presence of the vector field u"" in the 
structure. 

Here is an application of these ideas. Consider the geometrical structure 
consisting of a nowhere-vanishing vector field m", together with a nowhere- 
vanishing 3-form, uJabci that has zero curl and is annihilated by u"". This 
structure satisfies both of the features above — it is locally homogeneous, 
and it gives rise to a hyperbolization for (H), (pISl). So, this geometrical 
structure could serve as the basis for a Lagrange formulation. In fact, this 
formulation is appropriate for a physical system, namely that of a fluid with a 
2-dimensional manifold of internal states, as discussed in Sect. 3. Identify 
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with the velocity field of the fiuid, and uiabc with the particle-number density, 
via Uabc = neabcdu'^- 

It is curious that the original system and its Lagrange formulation, while 
so similar with regard to their solutions, are completely different with re- 
gard to their initial- value formulations. Indeed, as we have seen in Sect. 2, 
it is frequently the case that the original system of equations, (Q), has an 
initial- value formulation, while its Lagrange formulation, does not. 

Perhaps there is some more natural or more general notion of "initial-value 
formulation" that would resolve this disparity. 

Appendix A — Derivative Systems 

Fix, once and for all, a first-order, quasilinear system of partial differential 
equations, as described in Sect. 2. That is, fix a fibre bundle, with bundle 
manifold B, base manifold M, and projection mapping B ^ M, together with 
smooth fields k^°'a,j^ on the bundle manifold B. Our system of equations, 
on a cross-section, M B, of this fibre bundle, is given by 

fc^\(V0)„"=/. (30) 

We shall now construct from this system a new first-order, quasilinear 
system of partial differential equations. The idea is to "take one derivative" 
(with respect to the point of M) of Eqn. (^). 

The first step is to introduce the appropriate bundle of fields for the new 
system. Let the base manifold again be M. But now let the fibre, over a 
point X G M, consist of all pairs, {(pXa""), where is point of B satisfying 
7r(0) = X and (a" is a tensor at satisfying 

fc^\Ca"=/. (31) 

Thus, is merely a point of the fibre over x G M, in the original bundle 
B. It represents a set of "values for the original fields" at x. The tensor (a"' 
represents a set of "values for the derivatives of the original fields" . In order 
that a given (a°' be a viable candidate for these derivatives, it must satisfy 
Eqn. (|3T|), the algebraic equation that results from replacing (V0)a" in Eqn. 
([30| ) by (a"- We impose this algebraic condition on ( in the very construction 
of the new bundle (as opposed, e.g., to introducing it later as an "algebraic 
constraint"). In short, the dynamics (Eqn. (|30|)) of the original system goes 
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into the kinematics (Eqn. (0)) of the new system. Call the bundle space of 
this new fibre bundle B'. Thus, the dimension of the fibres of B' is given by: 
(dim fibres of B){1 + dim(M)) - (dim vector space of equations in B). 

Consider, as an example. Maxwell's equations. Then M is a 4-dimensional 
manifold, with fixed smooth metric Qab of Lorentz signature. For the bundle 
B, the fibre over x E M consists of all antisymmetric tensors, F^c, at x. Eqn. 



([30| ) is Maxwell's equations: g W aF^c = 0, V[aFbc] = 0. For this example, 
the new bundle, B', has, as its fibre over x G M, the collection of all pairs, 
{Fbc, Cabc), with symmetries Fbc = F[bc], (abc = (a[bc], and with ( satisfying the 
algebraic conditions (Eqn. (pi])) g°'^Cabc = 0,(iabc] = 0. Thus, the fibres of B 
have dimension six, those of B' dimension twenty-two. 

Returning to the general case, the second step is to introduce appropriate 
equations on this bundle. A cross-section of the bundle B' consists of fields 
0, Ca" on M. On such a cross-section, we impose the following system of 
partial differential equations: 

(V0)," = Ca", (32) 

V[„a]" = /afc". (33) 
Eqn. (p2D provides the "interpretation" of (, as the derivative of (p. The 



/ab" on the right of ( PB| ) is some field on B' (i.e., some function of (a;, 0, C)); 
whose exact form depends on what derivative operator is used on the left 
side of that equation. The general rule is that Eqn. (^) to be the result 
of taking the curl of Eqn. (0). For example, if is represented by tensor 
fields over M, if ( is represented by the tensor fields obtained by taking the 
covariant derivatives (with respect to some fixed derivative operator on M) 
of those fields, and if that same derivative operator is used on the left in Eqn. 



(p3|), then / will consist of certain terms involving and the curvature tensor 
of that derivative operator. If, on the other hand, all bundles are taken as 
simple products, and all derivatives are taken using the corresponding (flat) 
connection, then /ab" = 0. Note that we have not included in our system 
the derivative of Eqn. (PT|). The reason is that Eqn. ( PT| ) has already been 
included at the algebraic level in the construction of the bundle B'. Its 
derivative is thus an identity in B'. On the other hand, we do include in our 
system Eqn. (|33D, even though it merely results from taking a derivative of 
Eqn. (^2]). In short, all algebraic conditions on fields are included in the 



construction of the bundleQ, while all differential conditions on fields are 
^^In fact, there is, at this level of generality, a possible anomaly with the system 
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included in the equations on a cross-section of that bundle. We note that 
the system of Eqns. (|3^ ) -(|33| ) is indeed first-order and quasilinear. 

Consider again the example, above, of Maxwell's equations. Then a cross- 
section of bundle B' consists of smooth fields, Fbc, (abc, satisfying everywhere 
the symmetries and algebraic conditions given above. The equations, (|32D 
-(|33|), on such a cross-section become, respectively 



VaFbc = (abc, (34) 
V[rfCa]fec = 2-Rrfa[b'"Fc]m- (35) 

Given a system, consisting of bundle B and partial differential equations 
(|30|), then by its derivative system we mean the system, consisting of bundle 



B' and partial differential equations (|32D-(|33D, constructed above. Note that 
every solution of the original system gives rise to a solution of its derivative 
system (by merely setting (a"' = (V</')a"). Conversely, every solution of the 
derivative system gives rise to a solution of the original system (by merely 
ignoring (). The two systems of partial differential equations are, in this 
sense, "equivalent as to solutions" . But they are not "equivalent as to form" , 
a feature we exploit in Sect 2. 

We next turn to the issue of the existence of an initial-value formulation 
for these systems. As discussed in Appendix B, we say that a general first- 
order quasilinear system, (pO|), of partial differential equations admits an 
initial-value formulation provided it satisfies the following three conditions: 
i) the system admits a hyperbolization; ii) all constraints of the system are 
integrable; and iii) the system has the correct number of equations relative to 
the number of its unknowns. [See Appendix B for the details of what these 
conditions mean.] A key property of the derivative system is the following: 
// the original system, admits an initial-value formulation, then so does 
its derivative system, (^^j- j^ ). We check the three conditions in turn. 

Let the original system (jS^) admit an hyperbolization (say, hjSA, with Wa). 
Then, we claim, so does its derivative system. Indeed, the corresponding 
bilinear expression (on a pair of tangent vectors, represented as (50", 5Ca") 
and (5'0°) (^'Ca")) is given by 



'^'5Ca°5'C/ + SrS'4>% (36) 



331) . In some cases, further algebraic conditions on the fields can follow from Eqn. (pSl). 
In fact, this anomaly will never arise in systems of interest, because it is precluded by 
the requirement, which we shall impose shortly, that all constraints of the original system 
be integrable. 
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where 9 is any positive-definite metric field on M. It is apparently not 
known whether the converse is true, i.e., whether the existence of a hyper- 
bolization for the derivative system, (^)-(Q), implies existence of a hyper- 
bolization for the original system, (^). Simple examples suggest that this is 
a reasonable conjecture. 

Integrable constraints of the system ( pUD do not lead to constraints of 
the corresponding derivative system. Rather, they lead to a reduction in 
the number of effective equations. Indeed, let a be any constraint. Then 
the result of contracting Eqn. ( |33D with c^a^^^o is an identity: It holds 
automatically, by virtue of Eqn. (^T]). Thus, each constraint for the system 
(pop reduces by one the number of effective equations represented by Eqns. 

dH) -m- 

What, then, are the constraints of the derivative system (P^-(P5D? These 
fall into two classes. The first class consists of those constraints that corre- 
spond to taking the curl of Eqn. (^2]). These constraints are of course in- 
tegrable: Their integrability conditions are precisely (^3]). The second class 
of constraints consists of those that correspond to taking the curl of Eqn. 
(]33|). These constraints, too, are integrable, by virtue of the fact that Eqn. 
([33| ) is itself a curl. Not all of these constraints, it turns out, are in general 
algebraically independent. 

Let us return to our original partial differential equation, (|30D . Denote 
by n the dimension of the base manifold M (the "number of independent 
variables" ) , by m the dimension of the fibres in the bundle E (the "number of 
unknown functions"), and by e the dimension of the vector space of equations, 
([30|). Further, denote by c the dimension of the vector space of constraints, 
and, for fixed nonzero covector w^, by c the dimension of the space of vectors 
of the form Wad^A, as d^A runs over the constraints. Then, as discussed in 
Appendix B, the condition that the original system ( ^0]) have the "correct 
number of equations" becomes e — c = u. We turn now to the derivative 
system (|32D-(p3D. The number of its unknowns is given by m' = m + {nu — e) 
(the two terms representing the numbers of unknowns contained in the fields 
(f) and respectively). The number of its equations is given by e' = nu + 
[un{n—l)/2—c\ (the two terms representing the number of effective equations 
in (^) and (0), respectively). Finally, the number of effective constraints of 
the derivative system is given by c' = u{n — l) + [{n — l){n — 2)u/2 + c — c] (the 
two terms representing the number of effective constraints in (|32D and (^Sj) , 
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respectively^). From these formulae, it is easy to check: li e — c = u, then 
e' — c' = u' . In other words, if the original system has the correct number of 
equations, then so does the derivative system. 

We conclude, then, that, beginning with a system (^) having an initial- 
value formulation, its derivative system, (il)-(il), also has an initial- value 
formulation. 

The construction above of the derivative system is useful because it per- 
mits a large class of systems of partial differential equations to be cast into 
a form to which the Lagrange formulation of Sect. 2 can be applied. But, 
unfortunately, passing to the derivative system and then to the Lagrange 
formulation is often a cumbersome procedure. The reason is that the deriva- 
tive system requires the introduction of additional fields to represent the 
derivatives of all the fields of the original system — even of those fields only 
remotely related to the one real interest: the velocity field. The result is a 
large number of extraneous fields. More useful would be a construction that 
goes only part way to the full derivative system — one that introduces ad- 
ditional fields to represent the derivatives of only some of the original fields, 
leaving the remaining ones intact. It turns out that, while there are one or 
two systems (e.g., that for dust) for which a smaller derivative system along 
these lines is available, for the vast majority of systems of partial differential 
equations of physical interest there is none. Here, briefly, is why. 

First, we must designate which of the dependent variables (the fields 
represented by 0) are to be derived and which not. This is done by writing 
the original bundle, B, as a product of two bundles, B' and B" , with the same 
base spaced M. The bundle B' carries the fields whose derivatives will be 
represented by new variables, while B" carries the remaining fields. A cross- 
section (f) oi B consists precisely of a pair, (0', 0"), where 0' is a cross-section 
of the bundle B' , and (f)" is a cross-section of the bundle B" . In terms of these 
variables, Eqn. (|30|) becomes 

A;'^^'(V0')a"' + A:"^^"(V0")a"" = (37) 

^^The number of effective constraints of Eqn. (js^) is the dimension of the vector space 
of tensors K°-^a satisfying A'^''^ = Al'^^'l^ and Wah-^^a (namely, (n - l)(n - 2)u/2), 
minus the dimension of the vector space of such tensors of the form c°a^^''q for c"" a a 
constraint (namely, c — c) . 

^■^Recall that the product of two bundles, with the same base space Af, is the bundle, 
again with base space Af , whose fibre, over point x S Af, is given by the product of the 
fibres, over x, in the separate bundles. 
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where primed Greek indices denote tensors in B', and double-primed in B". 
Here, the fields k', k" and j are all functions on B, i.e., are functions of 
[x , (j)' , (j)") . We now proceed just as with the derivative system. Introduce a 
new fibre bundle, with base manifold again M, but with fibre over a; G M 
consisting of certain triples, (0', Ca" There must now be imposed on 

such triples all those algebraic conditions that fiow from (|37|). This is done 
as follows. At each point, denote by V the vector space of fiA satisfying 
HAk"^°'a" = 0. That is, V captures "those equations in (|37D that contain no 
derivative of 0"" . We now demand, in order that a triple (0', (a" , <P") give rise 
to a point of the fibre, the following: For every ha € V", fiAk'^°'a'Ca°' = I^A]"^- 
This is the fibre bundle for our new system. Let the equations of the new 
system be 

(V0)/ = Ca"', (38) 
V[„ar' = /afe"', (39) 
yAk'^\<a''' + ^/Afc"^^"(V0")a"" = ^Aj^- (40) 



In (^OD , ua is any vector in some fixed subspace complementary to the sub- 
space V . In other words, Eqn. (^) refiects those equations of ( PTj ) that do 
involve the derivative of (j)" . 

The system (p^)-(|^) is, certainly, a first-order, quasilinear system of 
partial differential equations; and it has as its variables precisely the ones we 
intended, namely (0', Ca" 5 0")- But, unfortunately, this system is subject to 
a variety of maladies — and these can arise even if the original system was 
quite well-behaved. For example — and this happens frequently — there 
can be constraints for the system (p8D-(^) that are hidden in Eqn. (^OD, 
and thus do not arise from any constraints for the original system, (^7]). 
Furthermore, these new constraints are not in general integrable. One could 
attempt to include the integrability conditions of these new constraints as 
new equations for the system. But two further problems can arise. First, 
some integrability conditions can turn out to be mere algebraic equations 
on the fields, (0', Ca" ,0")- The only way to "include" such equations is to 
start over, introducing a new bundle right from the beginning. Second, some 
integrability conditions can turn out to be quadratic, rather than linear, in 
the field-derivatives. These cannot simply be "included" — at least, not if we 
wish to retain a quasilinear system. The system (^)-(|40|) can also manifest a 
number of other types of difficulties, e.g., absence of a hyperbolization or the 
wrong number of equations. There appears to be no simple, general condition 
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that guarantees that Eqns. lead to a system with an initial- value 

formulation. 

As an example of this construction consider again the simple fluid, (|I])- 
(0). Let B' be the bundle whose fibre consists only of the variable u""; and B" 
the bundle whose fibre consists only of the variable p. In this example, the 
vector space V, capturing those equations in ([l|)-(0) involving no derivative 
of p, is zero-dimensional. The corresponding new bundle space, then, is that 
whose fibre, over x G M, consists of {u"" , Wb'^ , p) , with u"" unit timelike and 
Wb°' satisfying QacU'^Wh"' = (unit-ness of m"). The equations for the new 
system, in this example, are 

Vbu" = (41) 

V^aWb]' = RabdW, (42) 

(^"- + n'^M'") V^p + (p + p)u^w^'' = 0, (43) 
u^V^p+{p + p)w^'^ = Q. (44) 
This system has a new constraint (obtained by combining Eqns. (|^) and 



to obtain an expression for VmP, and then taking its curl), which turns 
out not to be integrable. But its integrability condition turns out to be 
quasilinear in field-derivatives, and so may be included as a further equation 
of the system. The resulting system in this case (but not for the case of an 
even slightly more complicated fluid) actually admits a hyperbolization. 



Appendix B — Initial- Value Formulation 

Consider a first-order, quasilinear system of partial differential equations, as 
described in Sect. 2. That is, we have a fibre bundle, with base manifold M, 
bundle manifold B, and projection mapping B —>■ M. The system of partial 
differential equations, on a cross- sect ion, M ^B, of this bundle, is given by 
Eqn. d^). We are concerned here with the issue of under what circumstances 
such a system admits an initial- value formulation, i.e., a formulation in which 
the fields are first specified on some "initial surface" in M, and are then 
determined elsewhere in M by Eqn. (^ itself. 

The key to achieving such a formulation is an object called a hyperboliza- 
tion of the system (|^), a field /i^^ on the bundle manifold B having the 
properties described below. Consider, for {x, 0) any point of the bundle 
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manifold B, Wm any covector at x G M, and 50", S'(f)" any two vectors at 
(x, 0) G i3 tangent to the fibre ("vertical"), the expression 

WmhAk^"'jrS'<p''. (45) 

We demand, in order that this h^A be a hyperbolization, that, everywhere in 
B, this expression be symmetric in 50", 5'0" for all tUm, and positive-definite 
(i.e., positive for any nonzero 5'0^ = 50^) for some Wm- The most direct way 
to specify a hyperbolization for a system of partial differential equations is 
simply to give the bilinear expression (^). Such an expression indeed defines 
a hyperbolization provided it is symmetric and positive-definite, as described 
above, and furthermore, that it is some multiple of the result of replacing, in 
the left side of Eqn. (^, "(V0)a"" by "Wa50"". As an example, consider the 
system, ([|)-(||), for a simple perfect fluid. Consider the bilinear expression 

6'u^[{p + p){u^Wm)gabSu'+idp/dp)wJp] (46) 

+ {dp/dp){p + pr'6'p[{p + p)6u'^Wm + u"'wmSp]. (47) 

We note that this expression is symmetric under interchange of the two 
vectors ((5p, 5m") and {6'p,6'u'^), and that (provided {p + p) > and 1 > 
{dp/ dp) > 0) it is positive-definite whenever Wm is future-directed timelike. 
Furthermore, this expression arises, as described above, from Eqns. (|I])-(0). 
This bilinear expression, then, specifies a hyperbolization for this system. 

Let there be given a hyperbolization, haA, for the system (^. Then this 
object gives rise to an initial-value formulation for a portion of that system, 
in the following manner. Fix initial data, consisting of a submanifold T of 
M of codimension one (an "initial surface") together with a cross-section 0^ 
over this submanifold ("data" on that surface), such that at each point of T, 
the normal to T is one of the vectors Wm for which the bilinear expression 
(|45| ) is positive-definite (the surface is "non-characteristic"). Then: In some 
neighborhood of the submanifold T, there exists one and only one solution 
of the system 

V^^^(V0)a" = (48) 

such that = 0o on T. Note that we do not guarantee a solution of the 
entire system (^, but rather only of those components that are involved in 
the hyperbolization. While the proof of this theorem is technically difficult, 
the key idea is to construct, using the hyperbolization, an energy integral, 
which is positive-definite, and, effectively, conserved. 
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Denote by u the number of unknowns of the system (j^) (i.e., the dimension 
of the fibres in B), and by e the number of equations (i.e., the dimension of 
the vector space in which the index "A" hes). Then the mere existence of a 
hyperbohzation for this system already imphes e > u (i.e., that there are at 
least as many equations as unknowns). Should it happen that this inequality 
is an equality, i.e., that e = u, then it follows that the hyperbohzation tensor 
haA is invertible, and so that the system (^) exhausts the original system 
of equations, (0). Thus, in this case we are done: We have achieved our full 
initial-value formulation. In the example of the simple perfect fluid above, 
for instance, we have e = u = 4, and so the hyperbohzation { ^7\ ) gives rise 
to an initial- value formulation for the fluid system (||)-(@). Unfortunately, 
in many cases of interest we have the strict inequality e > u, i.e., there 
are additional equations in (^ that are not accounted for in (|48|) . Such 
"additional equations" are dealt with in the following manner. 

By a constraint of the system, (|]), of partial differential equations, at a 
point of B, we mean a tensor c"^ at that point such that the tensor c'^A^^^a 
is antisymmetric in the indices ^^a,b" . This definition has two facets. First, 
each constraint gives rise to an integrability condition. Fix a constraint field, 
c°'A, and a solution of Eqn. (H). Contract both sides of Eqn. (|^) with c^^, 
and apply to both sides some derivative operator, Vb, on M. Then, by the 
constraint-condition, terms involving second derivatives of vanish, leaving 
an algebraic equation (indeed, a polynomial of degree at most two) in the 
first derivative, (V0)a", of 0. The constraint field is said to be integrable if 
this equation is an algebraic consequence of Eqn. (^, i.e., if the difference 
of its two sides is the product of some expression (at most linear in field- 
derivatives) and the difference of the two sides of (^. Lack of integrability 
of a constraint generally indicates that "not all the equations have been 
included in the original system (^". As to the second facet, each constraint 
gives rise to a compatibility condition on initial data. Fix constraint field, 
c°'A, solution of Eqn. (^, and submanifold T of M of codimension one. 
Then, at each point of T, we have 

n„c'"AA;^^(V0)„" = n^c^Aj'', (49) 

where Um is the normal to T at that point. But, by virtue of the constraint- 
condition, the index "a" in the tensor ^mC^A^^^^o is tangent to T. Thus, Eqn. 
(|49| ) takes the derivative of only in directions tangent to T, and so it refers 
only to the value of on T, i.e., only to the initial data on T. In short, Eqn. 
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represents a compatibility condition on initial data. If these compatibil- 
ity conditions were not satisfied, then we would have have no hope of finding 
a corresponding solution of Eqn. (|^). As an example, consider the Maxwell 
equation V[aFhc] = 0. This equation has a constraint. The corresponding 
integrability condition, obtained by taking the curl of this equation, is an 
identity, and so this constraint is integrable. The compatibility condition 
(|49| ) on initial data becomes, in this example, V ■ -B = 0. 

In the case in which e > u, i.e., in which the system (^) has more equations 
than unknowns, two further conditions must be imposed on the system. The 
first is that all the constraints be integrable. The second is that e — c = u, 
where c denotes the dimension of the vector space of vectors of the form 
WmC^A, for fixed Wm, as c"^a runs through all the constraints. This last 
condition means that any additional equations in @ that are not included 
already in (|^) are accounted for, effectively, by constraints. It states that 
d^) has the "correct number of equations" for its unknowns. In the case of 
Maxwell's equations, for example, all the constraints are integrable, as we 
have already remarked; and we have e = 8, c = 2, and u = 6, so there is 
indeed the correct number of equations. That is, the two further conditions 
above are satisfied in this example. 

Consider now a first-order, quasilinear system of partial differential equa- 
tions that satisfies the three conditions given above. That is, let the system 
i) admit a hyperbolization, ii) have all its constraints integrable, and iii) 
have the correct number of equations, as described above. It seems likely 
that such a system — possibly with some mild further conditions — must 
always manifest an initial-value formulation in some suitable sense. That is, 
we would expect that, given initial data for the system on a suitable sur- 
face T, satisfying on T the compatibility conditions (^9]), then there exists 
a unique corresponding solution of Eqn. (^) in a neighborhood of T. A key 
piece of evidence prompting this expectation is the following. There cer- 
tainly exists a solution of Eqn. ( ^HD manifesting the initial data, as we have 
already seen. Consider next the left sides of Eqn. (|49|) (as c"^a varies over 
all constraints). These expressions of course vanish on T, and, by virtue of 
the condition ii) and iii) above, satisfy a system of equations that express 
the "time- derivatives" (off T) of these expressions in terms of their "space- 
derivatives" (within T). Naively, we might expect that, as a consequence, 
these expressions must vanish in a neighborhood of T. But the vanishing of 
these expressions implies, again by condition iii) above, that Eqn. (|^) itself is 
satisfied everywhere in a neighborhood of T. Indeed, in all physical examples 
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of which we are aware — including all those discussed in this paper — this 
naive expectation is in fact borne out. Unfortunately, there is, apparently, 
no general theorem to this effect. Nevertheless, we shall, for convenience, use 
the expression "having an initial-value formulation" to describe systems of 
partial differential equations that satisfy the three conditions, i)-iii), above. 
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